Field enhancement in subnanometer metallic gaps 
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Motivated by recent experiments [Ward et a/., Nature Nanotech. 5, 732 (2010)], we present here 
a theoretical analysis of the optical response of sharp gold electrodes separated by a subnanometer 
gap. In particular, we have used classical finite difference time domain simulations to investigate 
the electric field distribution in these nanojunctions upon illumination. Our results show a strong 
confinement of the field within the gap region, resulting in a large enhancement compared to the 
incident field. Enhancement factors exceeding 10^ are found for interelectrode distances on the order 
of a few A, which are fully compatible with the experimental findings. Such huge enhancements 
originate from the coupling of the incident light to the evanescent field of hybrid plasmons involving 
charge density oscillations in both electrodes. 

PACS numbers: 78.67.Uh, 73.20.Mf, 73.21.Hb 



Introduction.- The study of the optical response of 
metallic nanostructures is revealing fascinating new 
physics. Special attention is being paid to the anal- 
ysis of the so-called optical and infrared gap nanoan- 
tennas, which consist of adjacent metallic segments, 
like nanorods, separated by a nanoscale gap.^^^ The 
ability of these systems to efficiently confine and en- 
hance optical fields is crucial for applications such as 
single-molecule surface-enhanced Raman spectroscopy 
(SERS)^^ or extreme-ultraviolet generations^ 

There are many techniques to measure the near- 
field distributions in gap nanoantennas such as scanning 
near-field optical microscopy^' and two-photon induced 
luminescence j^i^ However, these local probes have a lim- 
ited spatial resolution of 10 nm at best.^^ It would be thus 
highly desirable to develop new techniques or strategies 
that enable the extraction of information about the local 
fields in subnanometer metallic gaps, where the field en- 
hancements are expected to be largest. Very recently, 
we made a step in this direction, investigating the 
electronic transport through atomic-scale gold electrodes 
separated by a subnanometer gap under near-infrared 
laser irradiation, finding that the irradiation induces a 
d.c. photocurrent. By comparing this photocurrent with 
low- frequency conduction measurements, we were able to 
determine the optical voltage generated across the gap 
and, in turn, to infer the electric field in this region. En- 
hancement factors exceeding 10^ were reported, in line 
with previous estimates from surface-enhanced Raman 
measurements The goal of this Communication is 
to shed some light on the origin of this huge field en- 
hancement in these atomic-scale gap antennas. 

From a theoretical point of view, the field enhancement 
in finite metallic structures with nanometric gaps such 
as nanoparticle or nanorod dimers has been extensively 



studied (see Ref. [l3 and references therein). However, 
here we are interested in the analysis of a sub-nanometer 
gap formed between two atomic-scale electrodes, thus 
coupled to semi-infinite leads, where, to our knowledge, 
no related studies have been reported ^i^ii^ In this work 
we present an analysis of the optical response of gold 
atomic junctions with subnanometer gaps based on finite 
difference time domain (FDTD) simulations. Our main 
findings are: (i) field enhancements exceeding 10^ are 
possible for gaps of a few A,^^ which supports the main 
conclusion of Ref. 13; (ii) the huge enhancements are 
due to the excitation of a hybrid plasmon involving large 
localized charge distributions of opposite sign on either 
side of the junction, in analogy with the dipolar bonding 
dimer plasmons found in nanoparticle dimers; ^^^^^ and 
(iii) the plasmon resonances red-shift as the interelec- 
trode gap decreases, which also resembles the behavior 
found in nanoparticle^^ and nanorod dimers. 

System and methodology.- We model here the gold 
nanogaps of Ref. with the idealized geometry of 
Fig. [TJa). In this geometry, two extended gold tips end- 
ing in semispheres with a radius of 2 nm are separated 
by a distance d (this distance is assumed to be the actual 
distance between the electron clouds). We also consider 
that the junction is placed in vacuum and that there is 
no substrate. To analyze the optical response of this gold 
junction we have performed FDTD simulations where the 
main features to be considered are: (i) the two tip-like 
electrodes are coupled to infinite metallic surfaces, which 
are described by perfect metallic boundary conditions;^^ 
(ii) the structure is illuminated at normal incidence by a 
plane wave covering the whole region, and perfectly ab- 
sorptive boundary conditions are placed on the illumina- 
tion direction; (iii) the gold dielectric function used was 
extracted from ellipsometry measurements of a 20 nm 
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FIG. 1. (Color online) (a) Schematic representation of the gold junction considered in this work, (b) Spatial distribution of 
the field intensity (normalized by the incident one) in the x-z plane for an interelectrode separation d = 3 A. The wavelength 
of the illuminating light is A = 801 nm. (c) Field enhancement (E / Eq) evaluated in the middle of the gap as a function of the 
wavelength of the incident light for different interelectrode separations. In all cases the polarization is directed along the x axis 
(junction axis). The dashed and dashed-dotted lines are the corresponding results for a single tip at distances 0.5 A and 7.5 A 
from the tip along the x axis, respectively. The vertical dotted line indicates the wavelength 785 nm used in Ref. 13. (d) Field 
enhancement in the middle of the gap as a function of the interelectrode distance at resonance and at A = 785 nm. 



thin film;^^ and (iv) a non-uniform mesh with a smallest 
grid size of 0.25 A at the junction was used. All the sim- 
ulations were performed with the code FDTD Solutions 
(from Lumerical Solutions, Inc., Canada). 

Results and discussions. - In Fig.lTJb) we show the spa- 
tial distribution of the electric field intensity (normalized 
by the incident intensity) in the x-z plane of a junction 
with d = 3 A illuminated by an incident light of 801 nm 
with its polarization along the junction axis (x axis). The 
field is strongly localized in the gap region, and the inten- 
sity (i.e. \E\'^) is enhanced by a factor larger than 6 x 10^ 
in that region. We have systematically analyzed the field 
distribution as a function of both the interelectrode dis- 
tance d and the wavelength of the incident light. A sum- 
mary of the results can be seen in Fig. [U^c), where we 
show the field enhancement factor (E/Eq) in the middle 
of the gap on the junction axis, with Eq the amplitude 
of the incident field, as a function of the wavelength, and 
for different values of d ranging from 1 to 14 A. Notice 
the appearance of a resonance that shifts monotonically 
to the red as d decreases. At large separations, this reso- 
nance wavelength tends toward that exhibited by a single 
tip, around 600 nm (see dashed and dashed-dotted lines). 
On the other hand, the enhancement factor on resonance 
reaches values larger than 10^ for d < 3 A, which sup- 
ports the estimates reported in Ref. [l3|. 

The resonance seen in Fig. [TJc) clearly suggests that 
the incident light is exciting a plasmon-type mode. The 
comparison of the wavelength of that resonance with the 
one of a single tip and the fact the resonance red-shifts as 
d decreases both indicate that this mode can be consid- 
ered as a hybrid plasmon involving charge density oscilla- 
tions of opposite signs on both sides of the junction. The 
red shift is a simple consequence of the increasing inter- 



electrode interaction as d decreases, which in turn leads 
to a reinforcement of the electric field in the gap region. 
It is worth stressing that in this case the hybridization oc- 
curs between the continua of delocalized plasmon modes 
of the extended tips, rather than between localized plas- 
mons as in the case of subwavelength nanoparticles (see 
discussion below). 

The actual value of the field enhancement factor de- 
pends on the tip radius and it increases monotonically 
as the radius decreases. This is precisely the well-known 
"lightning-rod" effect, see e.g. Ref. '26'. In this case the 
lightning-rod effect cooperates with the plasmon excita- 
tion to greatly enhance the field locally at the region. 
To get an idea of the impact of this effect, we have re- 
peated the calculations with a smaller tip radius of 1.5 
A and found that, while the spectral response is almost 
identical, the field maximum increases by approximately 
a 16%. 

In Ref. 13 it was found that the enhancement decreases 
slowly with d (slower than 1/d). Our simulations show 
that this decay depends critically on the wavelength of 
the incident light. In Fig.lTJd) we show the enhancement 
factor in the middle of the gap as a function of d both for 
the resonant wavelength and for 785 nm, the wavelength 
used in Ref. 13. Notice that for the resonant condition, 
the field decays monotonically with the gap size approxi- 
mately as l/d^-^^, while for A = 785 nm the field does not 
decay for very short distances, and for d > 4 A decays 
slightly faster than 1/d'^. 

The results described above qualitatively resemble the 
predicted and observed optical response of nanoparticle 
dimers (see e.g. Ref. [2ll ) and nanorods dimers (see e.g. 
Ref. l23h . For instance, it is well-known that when two 
nanoparticles are placed next to each other, the plasmon 
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FIG. 2. (Color online) (a) Spatial distribution of the nor- 
malized field intensity in the gold junction for d = 6 A and 
A = 750 nm. (b) The same as in panel (a), but for two spheres 
of radius 2 nm separated 6 A and illuminated with light of 539 
nm. (c) Field enhancement (E / Eq) evaluated in the middle 
of the gap of the gold junction of panel (a) as a function of 
the light wavelength, (d) The same as in panel (c), but for 
the nanoparticle dimer of panel (b). 



modes of the individual nanoparticles interact, resulting 
in hybridized dimer plasmon modes whose energies can 
be strongly red-shifted with respect to that of the plas- 
mon modes of the individual nanoparticles. This is es- 
sentially what happens in the gold junctions considered 
here. One may then wonder to what extent the field en- 
hancement in gold junctions can be explained in terms of 
the physics of nanoparticle dimer s. To answer this ques- 
tion, we have considered a dimer comprising two gold 
spheres of radius 2 nm, which is the radius of the outer- 
most part of the gold tips.^^ In the upper panels of Fig. [21 
we compare the field distributions for a gold junction 
with d = 6 A (panel (a)) and a dimer with the same sep- 
aration (panel (b)). In each case the wavelength consid- 
ered is the one that gives rise to the maximum field (750 
nm for the junction and 539 nm for the dimer). Notice 
that while the symmetry of the field distribution in the 
gap region is practically identical in both cases, there is a 
big difference in the magnitude of the field enhancements 
(of more than an order of magnitude). Moreover, as we 
show in Fig. [2fc-d), there is also an important difference 
in the spectral response, where the resonance for the ex- 
tended electrode junction is considerably red-shifted as 
compared with the dimer case. The smaller field enhance- 
ment for the dimer is mainly due to its smaller scattering 
cross-section. Of course, in the dimer configuration one 
could reach the enhancement factors of the junction for 
the same gap size, but that would require the nanoparti- 
cles to have a considerably larger radius, which would 
be unrealistic for an atomic-size contact. Although less 
important, the difference in the field values reached at 
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FIG. 3. (Color online) (a) Spatial distribution of the normal- 
ized field intensity in the x-y plane for a gold junction with 
d = 6 A and illuminated with light of 495 nm with its polar- 
ization directed along the axis perpendicular to the junction 
axis (^/-axis). The wavelength used is the one at which the 
field reaches its maximum in the middle of the junction. 



the resonant conditions in both structures is also partly 
due to the frequency dependence of the gold permittivity, 
which in particular has a slightly larger imaginary part 
at the frequency of the plasmon resonance of the dimer. 

The plasmonic origin of the field enhancement suggests 
that it should be very sensitive to the polarization of the 
incident field. This is indeed the case, as we illustrate 
in Fig. [3l This figure shows the intensity distribution in a 
junction with electrode separation d = 6 A and A = 495 
nm, but this time the polarization is directed along the 
axis perpendicular to the junction axis (this figure has to 
be compared with Fig. [2fa)). As one can see, the field 
distribution is now quite different (with the charge oscil- 
lating back and forth in the transversal direction) and, 
in particular, the near-field is smaller than the incident 
field everywhere in the junction region. Notice that a 
field distribution like this does not generate an optical 
voltage across the junction and therefore, there would 
not be any photocurrent in a transport experiment like 
the one of Ref. 0. We have also verified that the strong 
polarization dependence of the field enhancement persists 
even in very asymmetric contacts with a pronounced mis- 
alignment of the axes of the tips. Further investigation 
into other tip geometries are ongoing. 

It is important to discuss the limitations of our clas- 
sical theory. First, we have assumed that the optical 
response of these nanocontacts is well described by a clas- 
sical frequency-dependent dielectric function. However, 
non-local effects can play some role in the outermost part 
of the electrodes, when the dimensions are smaller than 
the mean free path of the valence electrons that partic- 
ipate in the plasmon resonances^^ Second, the onset of 
quantum effects when the tips are very close can be much 
more important. As discussed in Ref. |32|, the onset of 
quantum tunneling and the concomitant charge transfer 
between the electrodes can lead to the appearance of a 
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charge transfer plasmon involving conduction electrons 
flowing back and forth between the metallic tips. A fl- 
nite electron density between the tips could give rise to 
a screening of the plasmonic interactions responsible for 
the strong red shift of the dipolar plasmons and, in turn, 
could reduce the fleld in the gap region^^ In any case, 
while our classical approach neglects detailed electronic 
structure, the essential physics of plasmon hybridization 
involving the continuum modes of the extended elec- 
trodes is expected to survive intact in realistic quantum 
mechanical calculations. 

Conclusions.- We have studied the optical response 
of atomic-scale gold junctions with subnano meter gaps, 
within the framework of classical electromagnetism. We 
have shown that the huge fleld enhancements reported 
experimentally^^ originate from the excitation of hybrid 
plasmons resonances involving charge oscillations in both 
electrodes. Such resonances red-shift as the gap size de- 
creases as a consequence of the increase of the interelec- 



trode interaction. Our results on the magnitude of the 
fleld enhancement clearly indicate that metallic nanogaps 
can be ideal templates for surface-enhanced Raman scat- 
tering of single molecules, very important for both molec- 
ular electronics and sensing applications.^^ 
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